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A new combinatorial, analytical approach to the ground- 
state energy problem of spin glasses with different concen- 
trations of ±J interactions is developed. The energy eo is 
expressed in terms of the fraction of broken bonds /lo and ex- 
panded into a fast converging series of the average length of 
a segment between two frustrated plaquettes of a minimum- 
weighted perfect matching A. The concept of so called dia- 
gram stability s is introduced in order to calculate coefficients 
of this expansion. Finally, the fraction as a function of the 
concentration p of negative bonds is obtained for triangular 
and square infinite lattices in excellent accordance with nu- 
merical simulations of large adequate systems. 



(1) 



i=0 



I. INTRODUCTION 

The Ising model on a regular lattice with ± J interac- 
tions between nearest neighbours is one of the simplest 
which is able to describe some physical features of spin 
glasses appearing in nature. Even for such "binary" in- 
teractions there is no analytical approach to the ground- 
state energy apart from unfrustrated Mattis spin glasses 
101 which can be transformed via gauge transformation 
into a ferromagnetic system. 

Exact, numerical methods aimed at finding the 
ground-state energy can only be applied to relatively 
small systems (as, e.g., the branch and bound method 
1^,^). Approximate methods are used for larger systems 
(as, e.g., matching algorithms However, the nu- 

merical effort increases at least polynomially with system 
size [|. Thus the precision of an extrapolation to the 
thermodynamic limit remains unknown. Therefore, any 
analytical or semi-analytical approach to the problem is 
worth the effort. The purpose of the present paper is to 
propose such an approach and to apply it to the square 
and triangular lattices in order to find the fraction of 
unsatisfied bonds in the ground state of an Ising system 
with randomly mixed ± J bonds. 

The immediate stimulation has come from the article 
of Liebmann and Schuster |^ who considered the aver- 
age minimal distance d{p) between a frustrated plaquette 
and the nearest other frustrated plaquette in the square 
lattice. They used the following expansion: 



where Ci is the conditional probability, namely the prob- 
ability that a plaquette is frustrated and all surrounding 
plaquettes up to the distance d = i are not frustrated di- 
vided by the concentration of frustrated plaquettes. The 
authors' main finding was that the parameter d{p) repre- 
sent a continuous function of the concentration p of nega- 
tive bonds. They stated that series (|l|) converged rapidly 
and found that an excellent approximation for d{p) was 
obtained using only the first three (or even two) terms. 
This allows one to hope that a similar approximation 
should also work for an analogous series expansion of the 
energy. 

The parameter d{p) was also investigated by Bendisch 
. He found that the maximum of d{p) at p « 0.10 
p 0.15) for the square (triangular) lattice reflects an 
inflection point in the gromid-state magnetisation. 

In the present work, a new method of combinatorial ap- 
proach to a minimum-weighted perfect matching of the 
graph of frustrated plaquettes has been developed in or- 
der to obtain the expansion of the ground-state energy 
as a function of the concentration of negative bonds. 



II. BASIC DEFINITIONS 

The Ising Hamiltonian without external magnetic field 
is given by: 



(2) 



The interaction between the i-th and j-th spins has a 
ferromagnetic character for Jy > and otherwise an an- 
tiferromagnetic one. We will use the term of a positive 
(negative) bond interchangeably with the term of a fer- 
romagnetic (antiferromagnetic) interaction. We restrict 
ourselves to the following probability distribution of Jij : 



/(Jy) = (1 -p)<5(Jy - 1) +p<5(Jy + 1) 



(3) 
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for nearest neighbours (otherwise f{Jij) = 0). The pa- 
rameter p G [0,1] has a clear interpretation as the con- 
centration of negative bonds or the probability that an 
interaction is antiferromagnetic. We will consider the two 
dimensional triangular and square lattices (coordination 
numbers 6 and 4, respectively). 

The unit cell of the 2D lattice will be called a plaque- 
tte. The term of frustration is reserved for systems: the 
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smallest system is a plaquette, so we will say that a pla- 
quette is frustrated (unfrustrated) when an odd (even) 
number of plaquette bonds is negative. This term is not 
adequate for bonds because there always exist two con- 
figurations of connected spins which locally minimize the 
energy. Thus, we will say that a bond is unsatisfied or 
broken when the configuration of connected spins does 
not minimize the energy locally. 

A useful formulation of the ground-state problem can 
be obtained using the duality concept |^ , ||] , as a prob- 
lem of finding a minimum-weighted perfect matching in 
the graph of frustrated plaquettes. Thus, the ground- 
state energy can be expressed as; 



Nc 



2wo 



(4) 



where N is the number of spins, c is the coordination 
number and Wq is the weight of a minimum- weighted per- 
fect matching which equals to the total number of broken 
bonds in the real lattice. The ground-state energy per 
bond Co reads: 



Co 



where the parameter 



fJ-o 



-1 + 2/io 



Wo 



Nc/2 



(5) 
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has a clear interpretation as the fraction of bonds un- 
satisfied in a ground state. The parameter /io admits 
an interesting generalization in jl^ where the so called 
"misfit parameter" is defined as: 



A=i 



where N\ is the number of segments of length A or 
Wo = ^A , 



(8) 
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with Nf as the number of frustrated plaquettes and A as 
the average segment length: 



A=^AP(A) 



(10) 
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P(A) is the probability that a randomly chosen segment 
has a length A. Clearly: 



P(A) = 



(11) 



and Nf/2 = E^A, so formula (|) holds. The misfit 
parameter can be written: 



'^°=A^^ 



(12) 



Elementary, geometrical considerations lead to a relation 
between the coordination number and the total number 
of plaquettes in a regular lattice (neglecting boundary 
effects): 
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(13) 
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min ailU ^jyidx 



are minimal and maximal energies, re- 
spectively, of an adequate "ideal" system which is anal- 
ogous to the original one but it is not frustrated. The 
advantage of definition (|^) is that it allows a compar- 
ison of frustrated systems of different nature. In our 
case, the "ideal" system is a Mattis spin glass, thus 
^min ~ ~^max ^-i^^ simple algebra leads again 

to formula (^. 

The notion of the fraction of bonds broken in a ground 
state appeared in the literature earlier jl^. For conve- 
nience justified by expression (||) for eo as a linear func- 
tion of /io, from now on we will be concerned with cal- 
culating the fraction /io, i.e., the misfit parameter. Also 
final results will be presented in terms of this quantity. 



III. SERIES EXPANSION OF THE MISFIT 
PARAMETER 

The weight of a perfect matching in a ground state 
is equal to the total sum of lengths of segments each of 
them connecting a pair of frustrated plaquettes: 



where Pf{p) is the probability that a plaquette is frus- 
trated in the system with a concentration p of negative 
bonds. From now on, we will omit in the notation Pf{p) 
the explicit dependence on p, and we will write just Pp. 
Finally, we obtain a useful formula: 
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At this point, we stress that the statistics P(A) de- 
pends on the choice of one of the ground states. For 
example, in Fig. |l| are presented two possible minimum- 
weighted perfect matchings: one of them provides two 
segments of length 2 and the other one gives one segment 
of length 1 and one segment of length 3. Thus, there is 
no universal distribution of P(A), but this distribution 
may vary from one ground state to another. The com- 
mon feature of these distributions is that they must give 
the same mean length of segments A. To keep symbols 
as simple as possible we will not introduce a new index 
like Pm(A), but the reader should keep in mind that cal- 
culations are performed for a specific ground state. The 
choice of this ground state will be discussed later. 
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IV. METHOD 



Referring to the discussion in section || we restrict our- 
selves to segments with lengths 1, 2 or 3. Hence, approx- 
imately, P(A > 4) = 0. It follows: 



^ ^ 3P(3) + 2P(2) 



P(l) 



P(3) + P(2) + P(1) 



A 



3,2-'l,2,l 



2P2 



1 



P3, 2^2,1 + R2,l + 1 



(15) 
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where P{i)/P{j). 

To perform calculations we choose a specific ground 
state, namely a state with as few as possible matching 
intersections. Fig. || illustrates the rule of thumb for such 
a choice (sometimes it is not possible to avoid intersec- 
tions - Fig. Ijc). Next we perform a transformation of the 
system such that in the transformed system only bro- 
ken bonds are negative (Fig. p|). This transformation is 



defined as the mapping: Ji^ 



and Ji. 



\Jii\ for broken bonds. 



I for satisfied bonds 
This mapping does 



not change the configuration of frustration of the system. 
Thus, both systems are equivalent with respect to the 
gauge symmetry in the meaning of [T^ . In particular, 
they have the same ground-state energy. Furthermore, 
in the transformed system the concentration of negative 
bonds is equal to the fraction /xq of bonds broken in the 
ground state, i.e., to the misfit parameter. 



A. i?2,i for the triangular lattice 

In Fig. ^, we show the three diagrams which contain 
segments of length 1 and the six diagrams with segments 
of length 2. In the present approximation, diagrams with 
matching intersections, such as shown in Fig. [^c are omit- 
ted. This is justified by our choice of the ground state 
with as few intersections as possible. Because of the ro- 
tational symmetry between diagrams considered one can 
write: 



D 6n2a „n2a 

^2,1 - t; = ^ 



(17) 



where nia and n2a are the probabilities of finding dia- 
grams of type la and 2a, respectively, in the lattice. 

To calculate the ratio n2a/nia we add a plaquette to 
the diagram la in order to have the same number of pla- 
quettes and bonds as in the diagram 2a. This additional 
plaquette may be either frustrated or not and we have 
three new diagrams: la', la", la'" (Fig. ||). If the lattice 
were a set of isolated diagrams (which is not the case) 
like diagrams la', la", la"' and 2a, we could write: 



»2a ^ P|(l-PF)A^g(l-M0)' 

n,a' P|(1-Pf)mo(1-A^o)^ 



Mo 



1 - A^o 



»2a 

nia" 

nia'" 



PF)yl{l - Mo)' 



PUI^Pf)h'o{^-1^o)' 



Pl{l-PF)^il{l 



Moj 



I-Pp 
Pp 

1 - Mo 
Mo 



(18) 



The equations above contain information on isolated di- 
agrams. Diagrams la' (or la'") and 2a have the same 
number of frustrated and unfrustrated plaquettes but dif- 
ferent numbers of satisfied and broken bonds. Diagrams 
la" and 2a have the same number of satisfied and broken 
bonds but differ in the numbers of frustrated and unfrus- 
trated plaquettes. All these diagrams have the same total 
number of plaquettes and bonds. Hence, for isolated di- 
agrams (i. e., independent of any infiuence of the outer 
part of the infinite lattice) (^ would hold. However, 
the diagrams are situated inside the lattice and feel the 
presence of other frustrations. Therefore, for diagrams 
within the lattice, we introduce a new parameter, called 
stability Sq of diagram a, such that: 



n2a 

nia' 

»2a 
nia" 

n2a 
nia'" 



M0S2a 



(1 - Mo)'Sla' 
(1 - PF)s2a 

PpSla" 
(1 ~ Mo)g2a 

MoSla'" 



(19) 



What is the meaning of the stability parameter? It 
measures the probability of the occurrence of such a frus- 
tration in the outer part of the lattice which destroys 
the matching of a diagram, i.e., some of the segments 
vanish and others appear. Fig. ^ shows an example of 
a destruction of the matching in diagram 2a. Destruc- 
tion of the matching appears if a loop exists containing a 
given diagram and for which an alternative matching is 
energetically more profitable. Accordingly, the stability 
parameter can be expressed as a loop series: 



^diag. 



E 

loops 



loop 



(20) 



where the sum goes over all loops containing a given di- 
agram and Dioop is the probability that there appears a 
frustration destroying the matching of a diagram inside 
a loop. 

Because a segment of length 1 is energetically the most 
profitable one, we admit: 



Sla' 



Sla" 



Sla" 



1 



(21) 



For S2a we will consider only the first order loops in the 
expansion (|2^). We are left with one type of smallest 
loops with six plaquettes (Fig. ||) . Frustration of plaque- 
ttes belonging to the diagram is fixed. Furthermore there 
are only two configurations of frustration of the remain- 
ing plaquettes which destroy the matching of diagram 
2a. One of these configurations contains two frustrated 
and one unfrustrated plaquette and the other contains 
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three frustrated plaquettes. Thus, the contribution to 
the stabiHty S2a up to first order is: 



and 



1 12 

De = -(1 - Pf)pI+pI = gPF + gPF 



Combining (|T|) and (|l|) we obtain: 



(22) 



(23) 



R2.1 = 2 



nia 



2 nia" 



nia'' 



n2a 



D. i?2,i for the square lattice 

As in the triangular lattice, we choose the ground 
state with as few matching intersections as possible and 
again we perform a gauge transformation such that in 
the transformed system only broken bonds are negative. 

There are two diagrams which contain segments of 
length 1 and six diagrams which contain segments of 
length 2 (Fig. 0). As before, we omit diagrams with 
matching intersections in accordance with the foregoing 
choice of the ground state. Thus, for the square lattice 
we obtain a formula analogous to (EtI): 



Pf 



R2 



Mo 



M0S2a (1 - PF)s2a (1 - ^^o)s2a 

_ 2^0(1 - /io)(l - PF)s2a 

" (1 - ^ioY{l - Pf) + 2/Xo(l - Mo)^F + ^2(1 - Pf) 

with S2a given by (ES 



(24) 



4n2Q 
2 nia 



2 n2e 

2 nia 



^n2a 



(28) 



B. ii3,2 for the triangular lattice 

There are six diagrams which contain segments of 
length 2 and twelve diagrams with segments of length 
3 (Fig. HI). As before, we omit diagrams with matching 
intersections. A practical approximation is made: 



R 



3,2 



12 n3a 
6 n2a 



'"la 



R2A 



(25) 



C. Final formula for the triangular lattice 

Using (|l|), (^ and (||) we get: 

3i?| 1 + 2i?2,i ' ^ 



2c^io - (c - 2)Pf 



-^2,1 



R 



2,1 



1 



= 



(26) 



where for the triangular lattice the probability that the 
plaquette is frustrated amounts to: 



Pf = p(4p2 - 6p + 3) 



(27) 



The last two equations combined with (g4|) for R2.1 and 
(^) for S2a make it possible to calculate the misfit param- 
eter ^0- To solve this system of equations numerically a 
simple bisection algorithm was applied. Results are pre- 
sented in Fig. [lC| a together with the data of simulations 
of large samples . 



However, not all segments of length 2 have the same sym- 
metry in the square lattice (Fig. ^ diagrams 2a and 2e), 
so unlike (|l^, eq. (|2^) is only an approximation. We add 
a plaquette to the diagram la to have the same number 
of plaquettes and bonds as in diagram 2a. We have three 
new diagrams: la', la" and la'" (Fig. ||). Applying the 
concept of the diagram stability for the square lattice 
we get an equation identical to (p^. As before, the ap- 
proximation (|l]) is made, since segments of length 1 are 
energetically the most profitable. 

Now, the first-order loop expansion ( po|) for S2a for the 
square lattice will be considered. We have two smallest 
loops with six plaquettes containing diagram 2a (Fig. |^). 
The contributions from these two loops are: 



Pf)p 



2 4 
+ 3^1^ 



Thus, for the square lattice we have: 
2 2 4 c. 

S2a ~ 1 - gPF - 3PF 



(29) 



(30) 



Formally, < Pp < 1 as Pp is a probability but for 
the square lattice Pp < 0.5, so there is no danger that 
thr above approximation leads to an unphysical, negative 
stability. Combining (|2|), (|l|) and (|j we obtain: 



i?2 



nia' 



3nia" + 3nia" 



n2a 



Mo 



Pf 



fJ-o 



AiOS2a 



(1 - PF)s2a 

3/io(l - /io)(l 



(1 - Mo)s2a 
PF)s2a 



(31) 



(1 - poYil - Pf) + 3po{l - /io)PF + 3/ig(l - Pp) ' 
with S2a given by (SO). 



E. i?3 2 for the square lattice 



There are six diagrams which contain segments of 
length 2 and fourteen diagrams with segments of length 3 
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(Fig. 0). We omit diagrams with matching intersections, 
as before, and make a practical approximation is made: 



-^3,2 ~ — 



14 7l3a 



6 n2a 



7n2a 7 
3 nia 9 



(32) 



F. Final formula for the square lattice 



Combining (|]J) and (g2|) we get: 

2i?2,l + 1 



7 jj2 

2cMo-(c-2)P;^^-^ 

9^2,1 



i?2 



1 



= 



(33) 



where the probabiHty that the plaquette is frustrated in 
the square lattice is: 



= 4p(l-p) (/ + (!_ p)2) 



(34) 



The last two equations when combined with (31) for R2.1 
and (|o|) for S2a make it possible to calculate the misfit 
parameter /iq numerically. Results are presented in Fig. 

® together with the data of simulations of large samples 
0- 



evident. We approximated the probability of finding a 
diagram or a loop by the formula (1 — Pp)^"p^' , where 
Nf is the number of frustrated plaquettes and is the 
number of those unfrustrated plaquettes out of which a 
given diagram is built. It is not quite correct and hope- 
fully could be replaced by an exact formula in future 
calculations. 

The method can be applied to any two dimensional 
regular lattice. An extension on the case of three dimen- 
sional cubic lattice is under the consideration. 
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V. CONCLUSIONS 



The ground-state energy eo per bond is a simple linear 
function (|^) of the misfit parameter /zq and the present 
approach provides values of the misfit parameter in the 
ground state of mixed Ising ferro-antiferro magnets in 
excellent agreement with numerical experiments for 
large triangular and square lattices. This proves the ef- 
fectiveness of the method and justifies a posteriori the 
approximations made. The results for the boundary val- 
ues of the concentration p of negative bonds, e.g., for 
p = or p = 1 are equal to the well known exact so- 
lutions in these cases. In fact, for the square lattice in 
both cases Pp = and because -R2.1 given by ( ^l| ) has 
no singularities the only solution of ( |33| ) is //q = 0- 

The situation for the triangular lattice is the same for 
p = 0. On the contrary, for p tending to unity Pp tends 
to 1 and ( |23| ) with ( p^ ) lead to i?2,i and /io — * 1/3. 
The plaquettes in the triangular lattice contain an odd 
number of bonds. Thus, in the antiferromagnetic limit 
the system becomes a fully frustrated one and elementary 
considerations prove that indeed one third of the bonds 
must be broken in any ground state (/uq = 1/3). 

The case of p = 0.5 is under particular considerations 
of many authors. Thus, the solutions for this concentra- 
tion of negative bonds are specified: p^^ — 0.2218 and 
Aio'' = 0.1546. 

Each step of the foregoing approximations can be de- 
veloped relatively easily up to the second and further 
orders of precision. The method is open to such improve- 
ments. At this point, one of possible directions of the de- 
velopment should be mentioned which is somewhat less 
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FIG. 1. Two different ground states contributing to differ- 
ent probability distributions P{A.). 
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FIG. 2. a) One specific ground state without a matching 
intersection of segments representing the basis of the method, 
b) A ground state with a matching intersection that can be 
transformed to a) via an inversion of the proper spin (marked 
with J,), c) A configuration of frustrated plaquettes necessi- 
tating intersection (disregarded in this work) . 

FIG. 3. The gauge transformation bound to a given ground 
state leading from the initial lattice (a) to the lattice with only 
broken bonds being negative (b). 

FIG. 4. Diagrams in the triangular lattice taken into ac- 
count during considerations. Diagrams with matching inter- 
sections were not considered. 

FIG. 5. Diagrams la', la", la'" and 2a taken into account 
in the calculation of the parameter R2,i of the triangular lat- 
tice. Note that there are two possibilities for diagram la" . 

FIG. 6. Two frustration configurations inside the smallest 
loop destroying the matching of diagram 2a of the triangular 
lattice. Solid lines represent new segments, crossed lines - 
destroyed segments. 

FIG. 7. Diagrams in the square lattice taken into account. 
Diagrams with matching intersections were not considered. 

FIG. 8. Diagrams la', la", la'" and 2a taken into account 
in the calculation of the parameter i?2,i of the square lattice. 
Note that there are three possibilities for diagram la" and 
also three possibilities for diagram la'". 



FIG. 9. The two smallest loops containing the diagram 2a 
in the square lattice which contribute to the stability calcu- 
lation. 



FIG. 10. Misfit parameter /ito as a function of the concen- 
tration p of negative bonds for the triangular lattice (a) and 
for the square lattice (b), obtained in this work. Crosses de- 
note results obtained by numerical simulations of large sam- 
ples (up to 200 spins) [13]. Rhombs denote average values 
over 500 samples for each concentration of negative bonds for 
square and triangular lattices of size 6x8 and 6x6, respec- 
tively [14]. 
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